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O ■ ABSTRACT 

<N ■ 

^ ■ We analytically and numerically study the motion of electrons along a mag- 

netic loop, to compare with the observation of the propagating feature of the 
CNJ ■ non-thermal microwave source in the 1999 August 28 solar flare reported by 

, , Yokoyama et al. (2002). We model the electron motion with the Fokker-Planck 

^ 1 equation and calculate the spatial distribution of the gyrosynchrotron radiation. 

^ ■ We find that the microwave propagating feature does not correspond to the mo- 

S-j ■ tion of electrons with a specific initial pitch angle. This apparent propagating 

feature is a consequence of the motion of an ensemble of electrons with different 
initial pitch angles, which have different time and position to produce strong 
radiation in the loop. We conclude that the non-thermal electrons in the 1999 
qq ■ August 28 flare were isotropically accelerated and then were injected into the 

loop. 

O ; 
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■ 1. Introduction 

X: 

In solar flares, a large amount of released energy is expended for production of non- 
thermal particles. They propagate in the solar corona and produce non-thermal emissions 
such as hard X-rays (HXRs), 7-rays, and microwaves. The observed non-thermal emissions 
have key information on acceleration and transport mechanisms of particles. 

Fl are non - therm al emissions are signatures of propagation of accelerated particles. For 



HXRs, ISakaol (119941 ) analyzed flares with double-footpoint HXR sources observed with 
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Yohkoh ( lOgawara et al.lll99ll ). and showed that the time variation of the emissions from 
double sources coincides with each other within 0.1 s. This indicates that electrons, not 
ions, propagate along the magnetic loop, reach the footpoints, and then emit HXRs there. 
Lin et al.l (120031 ) observed intense HXR and 7-ray emissions in the 2003 July 23 flare, b y 
using Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI; iLin et al.l 120021 ). 
They found that the location of the 7-ray line emission source, which corresponds to a pre- 
cipitation site of accelerated ions, does not coincide with the HXR sources. This implies that 
the acceleration site and/or propagation path of ions differs from that of electrons. 



For meter- wavelength radio emissions, the so-called type III radio burst (e.g-. lBastian et al 



19981 ) is thought to be due to electromagnetic conversion of Langmuir waves, which are ex- 
cited by electrons propagating from a flare site toward the interplanetary space along open 
magnetic field lines. For centimeter-wavelength radio emissions, the micro wave source lo 



cated at a remote distan ce from a flare main site (microwave remote source; lHanaokalll996 



19971 ; iNishio et al.lll997l ) is interpreted as the precipitation of electro ns propagating from the 
main site where the energy release and particle acceleration occur. lHanaokal (119991 ) showed 
that the light curve at the microwave remote source has a similar profile to the HXR light 
curve at the main si te, but is slightly delay ed by roughly the loop transit time of electrons. 
On the other hand, lAltyntsev et al.l (120081 ) reported a flare in which the microwave light 
cu rve slightl y leads the HXR. They argued that this time delay, opposite to that found 
by lHanaokal . could be realized if electrons are injec ted into an asymmetric magnetic loop 
toward a footpoint with a stronger magnetic field. iMinoshima et al.l (120081 ) analyzed the 
spatial and spectral variations of HXR and microwave emissions in the 2003 May 29 flare, 
and showed that the observed characteristics c an be interpreted as a c onsequence of the 
electron transport called trap-plus-precipitation (IMelrose fc Brown! 1 19761 ). 



There are a few microwa ye observations of the direct detection of propagating electrons 
from a time series of images. iBastian et al.l (119941 ) reporte d the motion of the 8.4 GHz and 
15 GHz microwave sources with a speed of 3,000 km s -1 . IWhite et al.l (120001 ) reported the 
rapid motion of the 0.33 GHz microwave source with a s p eed o f 2.6 x 10 4 km s _1 . The 



most striking observation was reported by lYokoyama et al.l (120021 ) (hereafter Y2002). They 
observed a flare that occurred on 1999 August 28, using Nobeyama Radioheliograph (NoRH; 
Nakajima et al.lll994l ) with 0.1 s cadence. From a time series of NoRH 17 GHz images, they 
found two different classes of propagating features of non-thermal microwave sources along 
the loop (see Fig. 3 in Y2002). One is a slower propagation at a speed of 6 x 10 3 km s" 1 , 
and the other is a faster propagation at a speed of 9 x 10 4 km s -1 . 



The origin of the slower propagating feature is discussed by IStepanov et al.l (120071 ). 
They interpreted this feature as propagation of low-frequency whistler waves that can trap 
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the parent electrons by wave-particle interactions at their turbulent front. The phase velocity 
of low-frequency whistler waves is roughly in agreement with the observation. The origin 
of the faster propagating feature is discussed by Y2002. They interpreted this propagating 
feature as a free stream of relativistic electrons along the loop. Assuming that the 17 GHz 
microwave-emitting electrons are almost relativistic (v ~ c > 9 x 10 4 km s _1 ), they concluded 
that the parent electrons are injected with large pitch angle, ~ 70°. 

The interpretation of Y2002 on the microwave propagating feature is subject to the fol- 
lowing assumptions. First, they assumed that electrons freely stream along the loop. How- 
ever, since the flare loop must be a converging magnetic loop, electrons suffer the magnetic 
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1981c |Pulkill985l ; lBastianl|l999l ). corresponds to the trajectory of the parent electrons. How- 
ever, since the gyrosynchrotron radiation mechanism intricately depends on many physical 
parameters, it is not necessarily evident that the apparent motion of the radiation source is 
identical to the trajectory of the electrons with a specific pitch angle. 

In this paper we reconsider the rapidly propagating feature of the microwave source 
in the 1999 August 28 flare and address the pitch-angle distribution of the electrons, by 
refining the treatments of the emission mechanism as well as the electron motion. In § [2] 
we present an analytic solution of the electron motion in a converging magnetic loop. § [3] 
shows our numerical model of the electron motion along the loop, in which we use the 
gyro-averaged Fokker-Planck equation to determine the electron phase space distribution. 
We calculate the gyrosynchrotron intensity distribution along the loop from the calculated 
electron distribution for comparison with the observation. In § H] we present our calculation 
results. The results do not support the interpretation of Y2002 but suggest that the electrons 
in the 1999 August 28 flare were isotropically accelerated and then were injected into the 
loop. In §[5] we summarize this paper. 



2. Analytic Treatment 

We show an analytic solution of the motion of a single electron in a converging magnetic 
loop that helps to understand calculation results in § [H This is described by the following 
equations: the equations of motion along the magnetic field line and the conservation of the 
magnetic moment, 



dr 

Tt = ^ 



(i) 
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B(r) 5(0) ' 1 J 

where r is the position measured from the loop top (r = is the initial position) along the 
field line, v is the velocity, /i is the pitch-angle cosine, Ho is the initial pitch-angle cosine, 
and B(r) is the magnetic field strength at the position, respectively. We do not consider the 
energy change of an electron, thus v = constant. 

For mathematical convenience, we give the magnetic field strength as a qua dratic func- 



tion with respect to r of the form (e.g., lAschwanden et al.lll996l ; lHanaokalll999l ): 



. 2 

r 




B(r) = B + (B t - B ) \-J =B < 

where B and Bf are the field strengths at the loop top and at the footpoint, rf is the 
footpoint distance from the loop top along the field line (a half length of the loop), and 
M = Bi/Bq is the magnetic mirror ratio, respectively. When its initial pitch-angle cosine 
Ho is smaller than the loss cone angle cosine fi c = a/1 — M _1 , an electron is trapped and 
has a bounce motion in the loop due to the magnetic mirror. With a help of equation (j3J), 
equations (OQ) and (121) give a solution of the bounce motion as simple harmonic oscillation: 

r = Rsm(ut), (4) 
H = Ho cos(atf), (5) 

where 

R = t =r f , (6) 

V(1-/^)(M-1) 

u = ^(I-hD(M-I)^. (7) 

Equation (JTj) tells that an electron with a smaller Ho (larger initial pitch angle) has a higher 
frequency of the bounce motion. 

We apply the above solution to the propagating feature of the microwave source reported 
by Y2002. We assume, in the same manner as Y2002, that the apparent propagating motion 
of the radiation source corresponds to the trajectory of electrons with a specific initial pitch 
angle. The observation gives constraints on several variables appearing in equations and 
(I7j). According to the analysis by Y2002 (see Fig. 3 in their paper), the electrons traveled 
~ 4.5 x 10 4 km from the injection site (loop top) to the loop end (footpoint). We adopt 
this value as the footpoint distance rf, and evaluate R ~ rf. An observed travel time of 0.5 
s is roughly equal to a quarter of the bounce period, u> ~ (27r)/(0.5 x 4) = it radian s _1 . 



- 5 - 



Assuming that the magnetic field strength at the emission site is on the order of 100 Gauss 
(by Y2002), th e velocity of non-thermal gyrosynchrotron-emitting electrons is almost the 
speed of light (IBastiaru Il999f ). «~c = 3x 10 5 km s _1 . Using these, we can evaluate the 
magnetic mirror ratio as well as the initial pitch-angle cosine from equations and ([7]), 

ujR , , 

/io = — -0.47, (8) 

M = 1 + ^4-1.3. (9) 

Y2002 obtained a magnetic field strength at the loop top, Bq ~ 200 Gauss by potential field 
extrapolation. The photospheric magnetic field strength at the loop end is ~ 400 Gauss. 
We consider that the electron propagation end is above the photosphere where the field 
strength is probably lower than that at the photosphere, because numerous collisions with 
ambient plasma prevent them from reaching the photosphere. Therefore we consider M < 2. 
From this inequality and equation ([7j), we obtain /i < 0.88. Consequently we evaluate 
0.47 < /i < 0.88, thus 28° < a < 62° which is smaller than that derived by Y2002. This 
is because we consider the field convergence while they did not. 



Numerical Model 



In the previous section we only treat the motion of a single electron. Our goal is to 
numerically model the gyrosynchrotron intensity distribution along the loop for comparison 
with the observation. For this purpose, we first numerically model the electron motion along 
the loop by the Fokker-Planck approach of the elect ron distribution in phase space. The 



gyro-averaged Fokker-P lanck equation is as follows (ILu fc Petrosianl Il988l ; lHamilton et al. 
1990l : lMcClementslll990h : 



dN n dN d lrt d 



EN 



dN 



(10) 



Here N(r, fi, E, t) is the electron distribution in phase space (number of electrons per unit 
length per unit pitch-angle cosine per unit energy), E and are the Coulomb energy loss 
rate and pitch- angle diffusion coefficient, /i is the magnetic mirroring force, E = T — 1 is 
the kinetic energy in units of the electron rest mass energy m e c 2 , V is the Lorentz factor, 
(3 = a/1 — r~ 2 , m c is the electron mass, respectively. 

We adopt the Coul omb energy loss rat e and pitch-angle diffusion coefficient for a fully 
ionized plasma given by iLeach fc Petrosianl (119811 ) , 



E = -Kn//3 = -u E E, (K = Aitcrl In A), 
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^ = ^(1-M 2 ) = p^(l-A (12) 



where r = 2.82 x 10~ 13 cm is the classical electron radius, n is the ambient plasma number 
density, In A is the Coulomb logarithm, and z/e = Kn/(f3E) is the Coulomb collision fre- 
quency, respectively. The ambient plasma number density and the Coulomb logarithm are 
treated as constant in this paper, n = 10 10 cm -3 and In A = 25. Since we address shorter 
timescale (order of 1 s) phenomena than the Coulomb collision time (order of 100 s for 1 
MeV electrons), the results and conclusion in this paper are almost independent of these 
variables. 

The magnetic mirroring force is given by 

1 „, ^dlnB(r) , , 

A = --c/3(l-/i 2 ) dr K) . (13) 

Knowledge of the field configuration is needed to give the strength Bir) in this equation, 
and also to allow the simulation of the spatial distribution of the emission, because the 
gyrosynchrotron radiation depends on a viewing angle with respect to the field line. We 
employ the two-dimen sional, stati c, and symmetric potential magnetic field in Cartesian 



coordinates as follows (jPriestlll982l ): 



B x = B { cos(kx)e kz , 
B z = — Bi sin(kx)e~ kz , 



(14) 



where x and z directions correspond to the tangential and normal directions relative to the 
solar surface, respectively. The parameter k is determined from rf and M (see below). This 
field configuration is shown in the left panel of Figure CD The strength of this field measured 
along a field line is written as 

B ^ = r- -if° h(h uv ( 15 ) 
cos [sin {tann^/crjjj 

where 

'1 + sin {cos" 1 (M^ 1 )} 



k — — In 
2r f 



sin {cos" 1 (M" 1 )} 



(16) 



The top right panel of Figure [T] shows the field strength measured along the field line. We 
use parameters rf = 4.5 x 10 4 km and M = 1.6 based on the estimation in § [2j The form of 
the field strength (eq. (|T5i) . solid line in the top right panel of Fig. CQ is very close to the 
quadratic function (eq. ([3]), dashed line), allowing us to use the analytic solution in §[5] as 
an approximation for the calculation result. 
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We give the initial condition for N as 
N(r, /I, E, 0) = exp 



0.2r f 



(17) 



where is the initial pitch-angle distribution given later. By this equation, energetic 
electrons with a power-law distribution are initially injected at the loop top. The initial 
condition is given to be symmetric in (r, /i) space, which yields the symmetric solution of 
equation ([10]): N(r,n,E,t) = N(-r, -ji, E,t). 

We numerically solve equation ffTOl) by using a finite difference method with operator 
splitting. The differential operators in equation ( jTOj) are split into four terms: the advection 
terms in r, /i, and E space, and the diffusion term. We solve the diffusion term by using the 
Crank-Nicholson method with central difference. We employ the symmetric boundary condi- 
tion in /i space to ensure total number conservation. For the ad yection term in E space we use 
an an alytic solution derived by the method of characteristics ( iCraig et al.lll985l : MacKinnon 
19861 ). and adopt a single power-law function for the necessary interpolation at the inter- 
mediate location between the grid points (if either grid point has a non-positive value, we 
instead apply a linear interpolati on) . For the advec tion terms in r and fi space, we adopt the 
CIP-CSL2 scheme developed by lYabe et al.l (1200 ll ). This scheme simultaneously solves the 
integrated value of N as well as N itself to keep "subgrid" information and to satisfy the mass 
conservation. A boundary condition at r = is d {N(r, fi, E,t) — N(r, —fi,E,t)} jdr = 0, 
for the symmetry. Another boundary condition at r = rf is as follows: (1) electrons with 
positive pitch angle, N(ri,fi > 0,E,t), are lost from the calculation box; and (2) there is 
no flow from outside the box, N(rt,n < 0,E,t) = 0. Physically, the former corresponds 
to the electron precipitation into the footpoint, and the latter means that the possibility of 
electrons originating below the footpoint (chromosphere) is excluded. We set 64 x 128 cells 
in (r, fi) space, and 128 logarithmically-spaced grids from 50 to 5000 keV in E space. 

Using the electron distribution N(r, fi, E, t), we further numerically calculate the distri- 
bution of the gyrosynchrotron radiation for comparison with the observation of Y2002. A 
general calculation of the gyrosynchrotron radiation in a magnetized plasma includ es effects 
such as self-absorption, absorption by ambient plasma, and the Razin suppression (IRamaty 
19691 ). These effects contribute at low harmonics {v/vb %j 10, where v-q is the electron gy- 
rofrequency) of the radiation in the optically-thick regime. However, Y2002 showed that the 
17 GHz emission observed with the NoRH in the 1999 August 28 flare is the optically-thin, 
non-thermal gyrosynchrotron radiation. Therefore, it is sufficient for our purpose in this pa- 
per to consider only the optically-thin radiation at high harmonics (1 < v/v-r% 100) from 
mildly relativistic electrons (r < 10). Under such limited condition, iPetrosianl (|198ll ) gave 
an useful approximation to calculate the gyrosynchrotron emissivity from mildly relativistic 
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electrons with arbitrary energy and pitch-angle distributions (eq. (5) in his paper). We 
adopt his formula to calculate the radiation from 17 to 34 GHz. 

The gyrosynchrotron emissivity j u (r,t) depends on a magnetic field strength and a 
viewing angle 9 with respect to the field line at the emission site as well as the parent electron 
distribution. We set a field strength at the loop top of Bq = 200 Gauss. A viewing angle, 
depending on the loop location, orientation, and tilt with respect to the Sun, is calculated 
from the direction of the field line and the line of sight. We consider two ideal cases of the 
viewing angle. When the line of sight is along the z-axis (see the left panel of Fig. [I]), 9 is 
given as a function of r, 

9(r) = cos -1 [- tanh(fcr)] , (18) 

that is shown in the bottom right panel of Figure [TJ We call this as the top- view case. Such 
situation is realized when the loop is at the disk center and is not tilted with respect to the 
solar surface. When the line of sight is along the y-axis, 9 = 90° is independent of r. We 
call this as the side-view case. Such situation is realized when the loop is on the limb with 
north-south orientation. 

We calculate the gyrosynchrotron intensity I u (r,t), which is the observed variable, from 
the emissivity for comparison with the observation. This can be obtained by considering the 
projection effect of the loop by an observer, 



L(r,t)=j u (r,t)^Y i (19) 
sm9 

where A(r) is the cross-sectional area of the loop determined from the conservation of the 
magnetic flux, B(r)A(r) = constant. 



4. Calculation Results and Interpretation 

We present our model calculation results and address whether the apparent propagating 
motion of the microwave source reported by Y2002 actually corresponds to the motion of 
electrons injected with a specific initial pitch angle. To discuss this, we perform calculations 
for two cases of the pitch-angle distribution of the initial condition (</>(/i) in eq. f|T7|) ): 
narrow-angle injection (§ 14.11) and broad-angle injection (§ I4.2p cases. 
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4.1. Narrow- Angle Injection 



We first consider that initially injected electrons have an almost unique pitch angle, 
following the interpretation made by Y2002. To simulate this case, we give an initial pitch- 
angle distribution peaking at /i = ±/x p : 



= exp 



0.1 



(20) 



The top left panel of Figure [2] shows the time variation of the electron number distribu- 
tion along the loop. We show the calculation result with fi p = y/l — M" 1 = 0.612, which is 
equal to the loss cone angle cosine /i c . The vertical axis corresponds to the spatial coordi- 
nates from the loop top to the footpoint. We use the 1 MeV electron distribution integrated 
over pitch angle for this illustration. The solid and dashed lines show the trajectories of 
electrons with initial pitch-angle cosine of respective ±0.61 and ±0.4, obtained by solving 
equation flTJ with equations (j2J) and (TToT) (not eq. ([3])) by using the 4th order Runge-Kutta 
method. 

Initially injected electrons move toward the footpoint and about a half of them reach 
there after 0.2 - 0.4 s. They are lost from the calculation box because they are initially in 
the loss cone. Remaining electrons bounce back by the magnetic mirror and move toward 
the opposite direction. Subsequent diffusion by Coulomb collisions is ignorable because the 
collision time (~ 200 s) is much longer than the loop transit time (< 1 s). 

We calculate the electron distribution not only in space but also in energy (see eqs. (TTOT) 
and (TTTj) ). Since the electron propagation is velocity-dispersive, the time variation of the 
electron number distribution differs with energy. However, relativistic electrons show almost 
the same profile as the top left panel of Figure [2] because their velocity is nearly c. 

The middle and bottom left panels of Figure [2] show the time variation of the 17 GHz 
gyrosynchrotron intensity distribution along the loop in the top- and side-view cases, respec- 
tively. Note that the vertical scale of the middle left panel is slightly different from other 
panels because of the projection effect of the loop. 

The strong radiation source is localized in space along the loop. In the side view (bottom 
left), radiation primarily comes from the footpoint. In the top view (middle left), strong 
radiation comes from the intermediate position between the loop top and footpoint. These 
results are attributed to the dependence of the gyrosynchrotron emissivity on a viewing angle 
as well as the parent electron distribution along the loop. We interpret these results by using 
the analytic solution in § [2J 
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The gyrosynchrotron radiation is primarily pro duced by mildly relativistic (~ 1 MeV ) 
electrons with the pitch-angle cosine of \x ~ (3 cos 9 (|Petrosianlll98ll ; iLu fc Petrosianlll990l ). 
Given all electrons with the initial pitch-angle cosine /io at the loop top, their pitch-angle 
cosine is given by equation ([2]) with as a function of r. The viewing angle 9 is also given 
as a function of r, which is different between the top- (eq. (TT8T) ) and side-view (6 = 90°) 
cases. By solving \i = (3 cos 9 with respect to r, we obtain the position r g at which strong 
radiation are produced. In the top view, the equation is written as 



f4 



[l-(4)(M- 1) (^j =/3tanh(> 



(21) 



The solution r g of equation (12 ip is smaller than rf, as long as /io is smaller than the critical 
value (obtained by taking r g = rf in this equation). In the narrow-angle injection case, we 
assume that injected electrons have an almost unique initial pitch angle, /io ~ /i p . When 
/i p = y/l — M -1 , equation (l2Tj) is 



(3 tanh (kr g 



\ 



M 



1- ^ 




(22) 



which gives r g of an intermediate value between and rf. Therefore the strong radiation 
comes from the intermediate position between the loop top and footpoint. This explains the 
intensity distribution in the middle left panel of Figure [21 



In the side view, 9 = 90° is independent of r. Then the equation /i — f3 cos 9 = results 



m 



/'o 



-rf, 



y/(l-ti)(M-l 

which again gives r g of an intermediate value between and rf as long as /io < /i 



(23) 



1 — M~ l . When /io ~ /i p = y 1 — M" 1 , this results in r g ~ rf. Therefore the strong 
radiation comes from the footpoint. This explains the intensity distribution in the bottom 
left panel of Figure [2j 

Based on these discussions, we conclude that electrons injected into the loop with an 
almost unique pitch angle do not yield the propagating feature of the radiation source along 
the loop. 



4.2. Broad-Angle Injection 



Next, we consider that initially injected electrons have an isotropic pitch-angle distri- 
bution: 0(/i) = constant. Electrons with a small initial pitch angle can reach the footpoint 
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while those with a large initial pitch angle are confined to a narrow region around the loop 
top. 

The top right panel of Figure [2] shows the time variation of the electron number dis- 
tribution along the loop. The solid, dashed, and dash-dotted lines show the trajectories of 
electrons with initial pitch-angle cosine of respective ±0.61, ±0.4, and ±0.2. As expected, 
electrons are broadly distributed in the loop compared with the narrow- angle injection case 
(top left panel). 

The middle and bottom right panels of Figure [2] show the time variation of the 17 
GHz gyrosynchrotron intensity distribution along the loop in the top- and side-view cases, 
respectively. The intensity distribution is quite different from that in the narrow-angle 
injection case (middle and bottom left panels). The strong radiation source is broadly 
distributed along the loop. We can see the propagating feature of the strong radiation source 
from the loop top to the footpoint in both the middle and bottom right panels, similar to the 
observation of Y2002. For example, rapidly propagating features are found during 0.3 - 0.5 
s and 0.9 - 1.2 s in the bottom right panel. These features do not result from the motion of 
electrons with a specific initial pitch angle but from the motion of an ensemble of electrons 
with different initial pitch angles. 

We discuss the apparent motion of the radiation source by using the analytic solution 
in § [2J For mathematical simplicity, we consider the side-view case. As mentioned in § 14. 1\ 
the gyrosynchrotron radiation is primarily produced by the electrons with \i ~ in the side 
view. It is found from equations (JSJ) and (J7J) that the pitch-angle cosine of an electron with 
the initial pitch-angle cosine /iq becomes zero at 



where m is an odd integer. This equation shows that when /i is smaller (i.e., an initial pitch 
angle is larger) t g is smaller. At t — t g , the electron position is given by equation (T4j) with 
sm(ujt g ) = 1. Using this, /io in equation (1241 can be removed and it is rewritten as 



This equation shows that when r g is smaller t g is smaller. Therefore, r g is smaller when (jlq is 
smaller. Electrons produce strong radiation at different time and position, depending on their 
initial pitch angle. By this equation, one can trace the position of the strong gyrosynchrotron- 
emitting electron in (r, t) space. Figure [3] is the zoomed image of the bottom right panel of 
Figure [2} to compare the analytic solution (eq. ( |25l) ) with the calculation result. The thick 



(m/2)vr 



(m/2)7T Tf 



(24) 




(25) 
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black line shows the analytic solution with m = 1 and v = c. As can be seen, the analytic 
solution well explains the propagating feature of the radiation source. 

Based on these discussions, we interpret the microwave propagating feature as follows. 
Electrons with a larger (smaller) initial pitch angle emit microwaves earlier (later) at the 
position closer to (farther from) the loop top, as seen in the peak of the dash-dotted (solid) 
line in Figure [H This difference, due to the difference of the initial pitch angle of the parent 
electrons, appears as the propagating motion of the strong microwave source along the loop 
from the loop top to the footpoint. 

The propagation of the radiation source seen in Figure [2] is periodic. This is because 
the initially injected electrons remain to oscillate in the loop (m = 1, 3, 5, . . .). Actual flares, 
on the other hand, show the light curve with a duration (tens of second) longer than the 
electron loop tran sit time scale ( < 1 s), indicating that electrons are continuously injected 



into the loop (e.g.. iMinoshima et al.ll2008l ). The periodic propagation seen in our calculation 
is probably occulted by the emissions from subsequently injected electrons, except the first 
one. The observable propagating motion is the first line (m = 1). The time of electrons 
reaching the footpoint and the propagation speed of the radiation source are calculated from 
equation ff25l) with r g = rf, 



7T / M r f n 2 M-l 

tcnd = 2 V M=T« ' Vpmp " W = -V ~ir V - (26) 

The propagation speed of the radiation source is determined from the electron velocity 
and the magnetic mirror ratio. When the magnetic field strength at the emission site is on 
the order of 100 Gaus s, effective en ergy of electrons producing > 10 GHz gyrosynchrotron 



radiation is > lMeV (jBastianlll999l ). hence their velocity is nearly c. The propagation speed, 
^prop — 0.4c with M = 1.6, is almost independent of the observation frequency, as long as 
parent electrons are relativistic. This frequency-independent microwave propagating feature 
is found not only at 17 GHz but also at higher frequencies up to 34 GHz in our calculation. 

We performed the calculation for two ideal cases of the view (top and side) because the 
gyrosynchrotron radiation depends on the viewing angle. The apparent propagating motion 
seen in the top-view case (middle right panel of Fig. [2]) looks different from that in the side- 
view case (bottom right panel), but can qualitatively be understood in the same manner 
as the side view. Difference of the view, which depends on the magnetic field configuration 
and location of the loop, does not influence our interpretation on the microwave propagating 
feature. 



It is supposed that a specific electron accounts for the propagating feature of the mi- 
crowave source. To produce the strong radiation, the pitch angle a of the electron should 
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be close to 9 at any time and position throughout its propagation along the loop. This is 
satisfied when (1) both a and 9 do not change, or (2) 6 changes in the same manner as a. 
The former case is realized if the field strength and the viewing angle do not change along the 
loop. The latter case may be realized if the loop is on the limb with east-west orientation, 
because 9 as well as a increase toward the footpoint in such configuration. These cases are, 
however, limited ones. Our interpretation holds in any probable field configuration and loop 
location. 



5. Summary and Discussion 

We presented analytic and numerical treatments of the electron motion along a magnetic 
loop, and consider the rapidly propagating feature of the non-thermal microwave source 
reported by Y2002. We studied this issue by describing the electron distribution with the 
Fokker-Planck equation and by calculating the gyrosynchrotron radiation. 

We first assumed that electrons injected into the loop have an almost unique initial 
pitch angle, following the interpretation made by Y2002. These electrons do not yield the 
propagating feature of the radiation source along the loop. This does not support the 
interpretation of Y2002. 

We next assumed that electrons injected into the loop have an isotropic pitch-angle 
distribution. In this case, the intensity distribution shows the apparent motion of the strong 
radiation source from the loop top to the footpoint, similar to the observation. This feature 
is interpreted as the motion of an ensemble of electrons (not a specific electron), which have 
different time and position to produce strong radiation due to the difference of their initial 
pitch angle. To show the microwave propagating feature, the injected electrons should be 
broadly distributed in pitch-angle space. 

We discuss a probable physical process of electrons in the 1999 August 28 flare. It is 
thought that this flare was triggered by the interaction bet ween the microw ave-propagating 



loop and a compact loop. S uch configuration is suggested by lHanaokal (Il999l ). called "double 



loop flare" . lHanaokal (119991 ) concluded that in the double-loop flare electrons are accelerated 
at the region where two loops interact. Based on this model, Y2002 interpreted that the 
acceleration site in the 1999 August 28 flare is where the propagation of the non-thermal 
microwave source starts, that is, the injection region. We conclude that the electrons in 
the 1999 August 28 flare started to propagate along the loop just after being isotropically 
accelerated at the site. 

In § |4.1l we discussed the calculation for the narrow-angle injection case with /i p = 0.612 
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only. The discussion can be applied to electrons with different pitch-angle distributions. 
When the electron pitch-angle distribution is concentrated perpendicular to the magnetic 
field, /i p = 0. In this case, it is found from equations (|2"T|) or (1231) with fi ~ that the strong 
radiation source is confined to the loop top. When the electron pitch-angle distribution is 
concentrated parallel to the field, on the other hand, fi p = 1. In this case the solution of 
equations (T2TT) or (1231) with fj, ~ 1 is r g = oo, meaning that there is no position within 
the loop to satisfy the condition for strong radiation. Neither of these cases results in the 
microwave propagation. 



We calculated the gyrosynchrotron intensity from the approximation of lPetrosianl (119811 ). 
which is valid when ambient plasma density at the emission site is so low that the absorptions 
and the Razin suppression are n ot important. We have to use a more exact formula of the 



gyrosynchrotron intensity (e.g., iRamatyl Il969l ; iFleishman Sz Melnikovl 120031 ) in case that 



microwaves are produced at the high-density site in which these effects are important. 

As well as the intensity, one can utilize the degree of circular polarization of the gy- 
rosyn chrotron radiation fo r diagnostics of the pitch - angle distribution of parent electrons 



[e.g., lAltyntsev et al.ll2008l ). IFleishman fc Melnikovl (120031 ) found from their numerical cal- 
culations that the degree of polarization increases when the electron pitch-angle distribution 
is more anisotropic and/or the line of sight is more parallel to rather than perpendicular to 
the magnetic field line. Propagating electrons in a converging magnetic loop have a more 
anisotropic distribution at the footpoint than the loop top. Furthermore, it is expected for 
the disk flares (top view) that the magnetic field line at the footpoint is close to parallel to 
the line of sight. Therefore, gyrosynchrotron emission is more likely to be polarized at the 
footpoint than the top of the electron-propagating loop. The 1999 August 28 flare certainly 
showed that the degree of polarization increases toward the footpoints from the loop top (see 
Fig. 2 in Y2002), agrees with the above statement. 

The propagating feature of the microwave source gives us great opportunities to study 
the gyrosynchrotron radiation mechanism and the electron transport, and to constrain the 
pitch-angle distribution of the injected electrons which is crucially important to understand 
the electron acceleration me chanism. As is evi dent from many hard X-ray observations show- 
ing footpoint sources (e.g., ISakaolll994l . Il999l ). propagation of electrons and the microwave 



source along the loop commonly occurs in the flare. However, detection of such phenomena 
is quite difficult. The 1999 August 28 flare is the unique event in that NoRH detected the 
microwave propagating featu re during its o bservational period until end of 2004 since its 
operation start in 1992 June (IShimojoll2007l ). This event was an extremely well resolved one 
in both space and time by NoRH. 



For further study, improvement of radio observatories is important. Observations with 
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high temporal ( < 0.1 s) and spatial (~ 1") resolutions more clearly resolve the microwave 
source because it propagates with a speed close to the speed of light along a loop with typical 
length < 100" ~ 7 x 10 4 km. Such observations should be implemented at frequencies 
greater than ~ 10 GHz which correspond to the optically-thin regime of the microwave 
emission in typical solar flares, to obtain the spectral property of non-thermal electrons. 
Since the gyrosynchrotron radiation depends on the viewing angle with respect to the loop, 
the center-to-limb variation of the microwave distribution is studied to addre ss the e l ectron 
dist ribution in the loop. S uch study has been carried out statistically by e.g.. iKosugil (119851 ) 
and lSilva fc Valentd (120021 ) with spatially-unresolved data. Statistical study on the center-to- 
limb variation of the spatial distribution of the microwave emission further gives constraints 
on the pitch-angle distribution of the injection electrons. Our numerical study will be of help 
to the future observational study for understanding the electron dynamics in solar flares. 
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Fig. 1. — Left: The potential magnetic field configuration described by equation (1T41) . The x- 
and z-axes correspond to the tangential and normal directions relative to the solar surface, 
respectively. The y-axis is perpendicular to the plane. We solve the Fokker-Planck equation 
(eq. (TTU]) ) along the thick line which has the half loop length rf = 4.5 x 10 4 km and the 
magnetic mirror ratio M = 1.6. Top right: The magnetic field strength measured along the 
field line (thick line in the left panel). The solid and dashed lines correspond to equations 
f fT5|) and ([3]), respectively. Bottom right: A viewing angle 6 with respect to the magnetic field 
line (thick line in the left panel) described by equation ffTBT) . The line of sight is assumed to 
be along the z-axis. When the line of sight is along the ?/-axis, 9 = 90° is independent of r. 
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Fig. 2. — Time variations of the electron number (top) and the gyrosynchrotron intensity 
(middle and bottom) distributions along the loop. The vertical axis corresponds to the spatial 
coordinates from the loop top to the footpoint. The left and right panels correspond to the 
results calculated for the narrow- and broad-angle injection cases, respectively. Top: The 
1 MeV electron number distribution along the loop. Relativistic electrons show almost the 
same distribution as that of 1 MeV electrons. Middle: The 17 GHz intensity distribution 
along the loop in the top view. Note that the vertical scale of these panels is slightly 
different from other panels because of the projection effect of the loop. Bottom: The 17 GHz 
intensity distribution along the loop in the side view. The solid, dashed, and dash-dotted 
lines show the trajectories of electrons with initial pitch-angle cosine of ±0.61, ±0.4, and 
±0.2, respectively. 
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Fig. 3. — Time variation of the 17 GHz gyrosynchrotron intensity distribution along the 
loop in the side view during 0-0.8 s, calculated for the broad-angle injection case. The solid, 
dashed, and dash-dotted lines show the trajectories of electrons with initial pitch-angle cosine 
of ±0.61, ±0.4, and ±0.2, respectively. The thick black line shows the analytic solution (eq. 

TO. 



